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In this paper, we extend the study of Drell-Yan processes with antiproton beams already presented 
in a previous work. We consider the fully polarized p^p 1 ' — ■> fi + /i~ X process, because this is the 
simplest scenario for extracting the transverse spin distribution of quarks, or transversity, which 
is the missing piece to complete the knowledge of the nucleon spin structure at leading twist. We 
perform Monte Carlo simulations for transversely polarized antiproton and proton beams colliding 
at a center-of-mass energy of interest for the future HESR at GSI. The goal is to possibly establish 
feasibility conditions for an unambiguous extraction of the transversity from data on double spin 
asymmetries. 

PACS numbers: 13.85.-t,13.85.Qk,13.88+e 



I. INTRODUCTION 

Building the nonperturbative structure of the nucleon bound state in QCD requires, first of all, the knowledge of 
the leading (spin) structure of the nucleon in terms of quarks and gluons. The transverse spin distribution (in jargon, 
transversity) constitutes the missing cornerstone that completes such knowledge together with the well known and 
• measured unpolarized and helicity distributions Q, 

From the technical point of view, the transversity is not diagonal in the parton helicity basis, hence the jargon of 
chiral-odd function. But in the transverse spin basis it is diagonal and it can be given the probabilistic interpretation 
of the mismatch between the numbers of partons with spin parallel or antiparallel to the transverse polarization of the 
Oh! parent hadron. In a nucleon (and, in general, for all hadrons with spin i), the gluon has no transversity because of 
the mismatch in the change of helicity units; hence, the evolution of transversity for quarks decouples from radiative 
gluons. But it also decouples from charge-even qq configurations of the Dirac sea, because it is odd also under charge 
conjugation transformations. In conclusion, the transversity should behave like a nonsinglet function, describing the 
distribution of a valence quark spoiled by any radiative contribution Q. The prediction of a weaker evolution of 
transversity with respect to the helicity distribution is counterintuitive and it represents a basic test of QCD in the 
nonperturbative domain. 

The first moment of transversity is related to the chiral-odd twist-2 tensor operator (7^75, which is not part of the 
hadron full angular momentum tensor Q. Therefore, the transversity is not related to some partonic fraction of the 
nucleon spin, but it opens the door to studies of chiral-odd QCD operators and, more generally, of the role of chiral 
symmetry breaking in the nucleon structure. In particular, another basic test of QCD should be possible, namely to 
verify the prediction that the nucleon tensor charge is much larger than its helicity, as it emerges from preliminary 
lattice studies 

From the experimental point of view, the transversity is quite an elusive object: being chiral-odd, it needs to 
be coupled to a chiral-odd partner inside the cross section. As such, it is systematically suppressed in inclusive 
Deep-Inelastic Scattering (DIS) [||. The first pioneering work about the strategy for its measurement suggested the 
production of Drell-Yan lepton pairs from the collision of two transversely polarized proton beams |?J . By flipping the 
spin of one of the two beams, it is possible to build a double spin asymmetry in the azimuthal angle of the final pair, 
that displays at leading twist the factorized product of the two transversities for the colliding quark- antiquark pair. 
This is the simplest possible scenario, since no other unknown functions are involved. But, in principle, the transverse 
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spin distribution of an antiquark in a transversely polarized proton cannot be large. Moreover, the combined effect 
of evolution and of the Soffer inequality seems to constrain the double spin asymmetry to very small values 0, . 

Alternatively, in semi-inclusive reactions the transversity can appear in the leading-twist part of the cross section 
together with a suitable chiral-odd fragmentation function [ld| . For 1-pion inclusive production, like in pp^ — * ttX or 
ep' — > e'irX reactions, the chiral-odd partner can be identified with the Collins function However, the situation 
is not so clear, since other competitive mechanisms (like, e.g., the Sivers effect |l2j ^ can produce the same single spin 
asymmetry when flipping the spin of the transversely polarized target. This happens because one crucial requirement 
is that the spin asymmetry must keep memory of the transverse momentum of the detected pion with respect to the 
jet axis, and, consequently, of the intrinsic transverse momentum of the parton. Several nonperturbative mechanisms 
can be advocated to relate the latter to the transverse polarization (see, among others, the Refs. 0, 0, 0, H^l. 
ultimately involving the orbital angular motion of partons inside hadrons 0, Il8t Il9| . 

Rapid developments are emerging in this field. In particular, new azimuthal asymmetries are being deviced to extract 
the transversity at leading twist while, at the same time, circumventing the problem of an explicit dependence u pon 
the transverse momentum. When the 2-pion inclusi ve p roduction pfij is considered both in hadronic collisions |2l| 
and lepton DIS with a transversely polarized target |22l l2^ | , the chiral-odd partner of transversity is represented by 
the interference fragmentation function [24[, of which a specific momentum enters the leading-twist single spin 
asymmetry depending only upon the total momentum and invariant mass of the pion p air |25l |26| . 

Experimentally, some recent measurements of semi-inclusive reactions with hadronic |27| and leptonic |28| beams 
have been performed using pure transversely polarized proton targets. New experiments are planned in several 
laboratories (HERMES at DESY, CLAS at TJNAF, COMPASS at CERN, RHIC at BNL). In particular, we mention 
the new project of an antiproton factory at GSI in the socalled High Energy Storage Ring (HESR). In fact, the option 
of having collisions of (transversely polarized) proton and antiproton beams should make it possible to study single 
and double spin asymmetries in Drell-Yan processes [2^, [3(], EI HE 01 with the further advantage of involving 
unsuppressed distributions of valence partons, like the transversely polarized antiquark in a transversely polarized 
antiproton. 

In a previous paper |35|. we have explored the Drell-Yan processes ppW — ► /i + /i~ X. For the single-polarized one, 
in the leading-twist single spin asymmetry the transversity happens convoluted with another chiral-odd function |36| | , 
which is likely to be responsible for the well known (and yet unexplained) violation of the Lam- Tung sum rule, an 
anomalous azimuthal asymmetry in the corresponding unpolarized cross section |37l l38l l39j . Monte Carlo simulations 
have been performed for several kinematic configurations of interest for HESR at GSI, in order to estimate the 
minimum number of events needed to unambiguously extract the above chiral-odd distributions from a combined 
analysis of the two asymmetries. 

In this paper, we will extend that work by considering numerical simulations for the fully polarized Drell-Yan 
process p^ p^ — > /i + n~ X, again at several kinematic configurations of interest for the HESR project at GSI. Since 
this is the simplest possible combination at leading twist involving, moreover, the dominant valence contribution of 
the transversity, the goal is to possibly establish feasibility conditions for its unambiguous direct extraction from data 
on double spin asymmetries. Most of the details of the simulation have been presented in Ref. [3{| and they will be 
briefly reviewed in Sec. ^ together with a description of the kinematics. Results are discussed in Sec. IHII and some 
final conclusions are drawn in Sec. IIVI 




FIG. 1: The Collins-Soper frame. 



II. THEORETICAL FRAMEWORK AND NUMERICAL SIMULATIONS 

In a Drell-Yan process, a lepton with momentum k\ and an antilepton with momentum ki (with kf/ 2 ) ~ 0) are 
produced from the collision of two hadrons with momentum Pi, mass Mi, spin Si, and P 2 , Mi, S%, respectively (with 
P\{2) ~ ^1(2)1 ^1(2) = — 1' -^1(2) ' ^1(2) = 0)- The center-of-mass (cm) square energy available is s = (Pi + P 2 ) 2 
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and the invariant mass of the final lepton pair is given by the time-like momentum transfer q 2 = M 2 = {ki + k 2 ) 2 - 
In the kinematical regime where M 2 ,s — ► oo, while keeping the ratio < r = M 2 / s < 1 limited, the lepton pair 
can be assumed to be produced from the elementary annihilation of a parton and an antiparton with momenta p\ 
and P2, respectively. If and P 2 _ are the dominant light-cone components of hadron momenta in this regime, 
then the partons are approximately collinear with the parent hadrons and carry the light-cone momentum fractions 
< x\ — pf J Pi , X2 = P2 I 'P 2 ^ 1j with q + = pi, q~ = p^ by momentum conservation [3rj . It is usually convenient 
to study the problem in the socalled Collins-Soper frame [4(J] (see Fig. where 

i = 
z = 



and q r is the transverse momentum of the final lepton pair detected in the solid angle (0, (f>). Azimuthal angles are 
measured in a plane perpendicular to z,t, and containing h. 
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FIG. 2: The leading-twist contribution to the Drell-Yan process. 
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A. Double spin asymmetry 



If the invariant mass M is not close to the values of known vector resonances, under the above mentioned conditions 
for factorization the elementary annihilation can be assumed to proceed through a virtual photon converting into the 
final lepton pair. Then, the leading-order contribution is represented in Fig. |21|3- The hadronic tensor is given by 



W^=IJ dpidp+ Tr P u Si) ^ $ (pa; P 2 , S 2 ) l v ] 



q <-> -q 
fj, <-> v 



(2) 



where the nonlocal correlators for the annihilating antiparton (labeled "1") and parton (labeled "2") are defined as 



#(pi;i\,Si) 

HP2-,P2,S 2 ) 



(2tt)4 
d 4 z 



(PiSi\i>(z)TP(0)\Pi,Si) , 



(27T) 



je l ^ z (P 2 ,S 2 \ip(0)ip(z)\P 2 ,S 2 



(3) 



They correspond to the blobs in Fig. [5] and contain all the soft mechanisms building up the distribution of the two 
annihilating partons inside the corresponding hadrons. 

Inserting into Eq. J3J) the leading-twist parametrization for $ and $ in terms of the (un)polarized partonic distribu- 
tion functions |4l| . the Drell-Yan differential cross section for transversely polarized hadrons, after integrating upon 
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dq T , becomes [42 



dx\ dx 2 d£l 12q 2 



l + C osH)^2e 2 f fl( Xl )f{(x 2 ) 
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+ sin 2 9 cos 2(j) 

+|S T1 ||S„| sin 2 9 cos(20 - Si -0 S J e}h{{x 1 )h{(x 2 ) 

f 



+ (1 «-» 2) 



(4) 



where a is the fine structure constant, dfl = sin 9ddd<f>, e/ is the charge of the parton with flavor /, and S is the 
azimuthal angle of the transverse spin of hadron i as it is measured with respect to the lepton plane in a plane 
perpendicular to z and i (see Fig. The function // (x) is the usual distribution of unpolarized partons with 
flavor /, carrying a fraction x of the unpolarized parent hadron; h{(x) is the transversity for the same flavor and 
momentum fraction (analogously for the antiparton distributions). The function v contains the contribution of the 
parton distribution 36], which describes the influence of the (anti)parton transverse polarization on its momentum 
distribution inside an unpolarized parent hadron: it is believed to be responsible for the observed anomalous azimuthal 
asymmetry in the unpolarized Drell-Yan cross section, the socalled violation of the Lam- Tung sum rule |37l l38t l39| , 
which no QCD calculation is presently able to justify in a consistent way [43l l44l l4 q . 
The double spin asymmetry is defined as 



-4, 



da TT - da n 



sin 



1 + cos 2 



cos(20 - <j> Si - (f> S2 ) 



J2 f e)h{{ Xl )h{{x 2 ) + {1^2) 
E/ e}f{( Xl )f{(x 2 ) + (l<->2) 



(5) 



In the next Section, we describe the details for numerically simulating A TT using the cross section Q. 



B. The Monte Carlo simulation 



In this Section, we discuss numerical simulations of the double spin asymmetry A TT of Eq. (JSJ for the p^p^ — ► /i + fj,~ X 
process. Our goal is to explore if it is possible to establish precise conditions in order to determine the feasibility 
of an unambiguous extraction of the transversity from data. Most of the technical details of the present simulation 
are mutuated from a previous work, where we performed a similar analysis for the unpolarized and single-polarized 
Drell-Yan process. Therefore, we will heavily refer to Ref. |3f| and references therein, in the following. 

As for the kinematics, several options have been considered in Ref. 35j|. Since the cross section decreases for 
increasing t, events statistically tend to accumulate in the phase space part corresponding to small r, where x p (xp) 
fall into the range 0.1-0.3 dominated by the valence contributions. In fact, the phase space for large r = x p xp is 
scarcely populated because the virtual photon introduces a 1/M 2 oc 1/r factor and the parton distributions become 
negligible for x p { X p) — * 1. Since the spectrum of very low invariant masses contains many vector resonances, where 
the elementary annihilation cannot be simply described by a diagram like the one in Fig. [21 it seems more convenient 
to reach such low values of r by adequately increasing the cm square energy s [35l |. Moreover, in this case the 
elementary annihilation should not be significantly affected by higher-order corrections like subleading twists, and the 
leading-twist theoretical framework depicted in Sec. Ill Al should be reliable. Hence, for the HESR at GSI the most 
convenient setup seems to be the option where antiprotons with energy Ep collide against protons with energy E p . 
Nevertheless, as in Ref. |35| we will explore also the option where the antiproton beam with the same energy hits a 
fixed proton target. 

As for the collider mode, neglecting hadron masses we have 

s = {Pp + Ppf w 4E P E P , (6) 

because for the two colliding beams P p = —P p . As in Ref. [35|. we will select the kinematics where antiprotons have 
energy E p = 15 GeV and protons E p = 3.3 GeV, such that s ~ 200 GeV 2 . If M is constrained in the "safe" range 
4-9 GeV between the cc threshold and the first resonance of the T family, then r falls into the statistically significant 
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range 0.08-0.4 where the parton distribution functions are dominated by the valence contribution. At the same cm 
energy, we will consider also the range 1.5 < M < 2.5 GeV between the <f> and J/ip resonances, which corresponds to 
the even lower range 0.01 < r < 0.03. 

For the fixed target mode, which at present is the selected setup by the PANDA collaboration at HESR at GSI |46[ . 
we have approximately 

s= (P P + Pp) 2 ~2M p E p , (7) 

such that for the considered Ep = 15 GeV it results s w 30 GeV 2 . For this case, we will restrict the invariant mass to 
the range 1.5 < M < 2.5 GeV, corresponding to 0.07 < r < 0.2. In fact, for M > 4 GeV most events are characterized 
by large partonic fractional momenta, where the quark-antiquark fusion model cannot explain the dimuon production. 
The Monte Carlo events have been generated by the following cross section [35| : 



da „ 1 
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il[)jJr iIt(I K-A(q T ,x F ,M)F(x F ,T)2^c t (q T ,x F ,T)S l (e,(t),4> Sp ,4> Sp ), (8) 



or, equivalcntly, 



= K-(xp + x p )A(q T ,Xp-x p ,M)F'(xp,x p ) V c-(q T , Xp, x p ) 0, (j) (j) ) , (9) 
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where the invariant x F = x p — x p is the fraction of the available total longitudinal momentum carried by the two 
annihilating partons in the collision cm frame. As it has been stressed in Ref. |35j ]. the range of values for x F depends 
on the energy. For the considered case of E p = 15 GeV, it results —0.9 < x F < 0.9, where positive values correspond to 
small 9 angles in the Collins-Soper frame (see Fig.QJ, and viceversa. Equations JSJ and © imply the approximation of 
a factorized transverse-momentum dependence, which has been achieved by assuming the following phenomenological 
parametrization 

A(q r ,x„M}= >• L ; j „ , (10) 
Mf)l 

where a(x F ,M), b(x F ,M), are parametric polynomials given in Appendix A of Ref. 39] and q T — |q T |. Actually, the 
Drell-Yan events studied in Ref. |3|j were produced by it — p collisions; however, the same analysis, repeated for p — p 
collisions [i^. gives a similar distribution for q T not very close to and not much larger than 3 GeV/c. In addition, 
the above q T distribution was fitted for M > 4 GeV and is singular near M — 1.5 GeV j3^; for M < 4 GeV we have 
assumed the same form with the mass parameter in the coefficients a and b fixed to the value M = 4 GeV. 

In order to simulate Eq. Q, the events produced by Eqs. JSJ or © have been integrated in dq T . However, apart 
from theoretical problems related to unwanted soft mechanisms, it is anyway not possible to collect events with very 
small (/^ because of the collider configuration. Hence, events have been selected with q T > 1 GeV/c. In our previous 
work 35], the searched asymmetry was emphasized by this threshold. Here, no such enhancement can be present, 
of course, because of the further integration. But the unavoidable dependence in q T , introduced by the lower cutoff, 
reflects in a drastic variation of the size of the sample. A precise answer depends on the experimental setup, but for 
q T > 1 GeV/c approximately 50% of the initial sample is excluded, while for q T > 0.5 GeV/c this fraction of events 
is reduced to 20%. 

The experimental observation that Drell-Yan pairs are usually distributed with q T > 1 GeV/c |39j, suggests that 
soft mechanisms are suppressed, because confinement induces much smaller quark intrinsic transverse momenta, but 
for the same reason it also indicates sizeable QCD corrections to the simple parton model. QCD corrections in the 
Leading-Log Approximation (LLA) 0] would imply a logarithmic dependence on the scale M 2 inside the various 
parameters entering the parton distributions 0] contained in Eqs. JSJl and such that it would determine their 
DGLAP evolution. However, it must be stressed that the key scale is M, and its range here explored is the same of 
Refs. [3^,23 1 where the functions F and Cj in Eq. © [or F' and c' ; in Eq. @] are assumed to be independent of M, 



F'(x p ,x p ) = —Y^ e ) f{(xp) f{(xp) + (p~p), (11) 



In particular, F' is given by 
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i.e. it is the azimuthally symmetric unpolarized part of Eq. @ which has been factorized out. The unpolarized 
distribution // (x) for various flavors f = u,d,s, is parametrized as in Ref. ^t|. QCD corrections in the Next-to- 
Leading-Log Approximation fNLLA) 14a are responsible for the well known K factor, which is roughly independent 
of x F and M 2 but it grows like yfr |39l |. In accordance with Ref. j3f|, for the range of interest 0.08 < r < 0.4 we 
assume as the best compromise the constant value K = 2.5. But we observe that in an azimuthal asymmetry the 
corrections to the cross sections in the numerator and in the denominator should compensate each other; indeed, the 
smooth dependence of the spin asymmetry on NLLA corrections has been confirmed for fully polarized Drell-Yan 
processes at high cm square energies 

The whole solid angle (9, </>) of the final muon pair in the Collins-Soper frame is randomly distributed in each 
variable. From Eq. the explicit form of the q T -integrated angular distribution is 

4 

^2 c'^xp^p) S l (9,(/),(f) Sp ,(j) Sp ) = l + cos 2 



i=l 

+ v(Xp^Xp)_ sin 2 9cQs2(f) 

+ \ S t p I |S Tj , I c 4 (x p , Xp) sin 2 6 cos(20 - <j) Sp - <f> s _ ) . (12) 

Recalling that the azimuthally symmetric unpolarized part F'(xp, x p ) of the cross section has been factorized out, the 
functions v and C4 turn out to be 



v{Xp, Xp) 



v(Xp, Xp) 



a(xp,xp) - — - f - ■ - -. (13) 



E/ e ) Ii( x p) fl ( x p) + (p^p) 

E/ e / h{{xp)h{(x p ) + (p <->p) 

E/ e / I{( x p) f{( x p) + (p^p) 



As for the former, the corresponding azimuthal asymmetry has been studied in Ref. 1351 adopting the simple 
parametrization of Ref. [3fil | and testing it against the previous measurement of Ref. |39j|. The latter has been 
further simplified by assuming that the contribution of each flavor / to the parton distributions can be approximated 
by a corresponding average function [3{| : 

r J r T \.._. (Msg)) fol(sp)) ml 
{jl( x p)) \M X P)) 

Four types of analytic dependences will be explored for the ratio (h±(x)) / (fi(x)) , namely the constants 1 and 0, and the 
ascending and descending functions y/x and y/1 — x, respectively. All the functional forms satisfy the Soffer bound 
across the whole range < x < 1, and they show different behaviours in the most relevant range 0.1 < x < 0.4 for 
the considered kinematics, as it will be clear in the following. The goal is to explore under which conditions such 
different behaviours in the ratio can be recognized also in the corresponding double spin asymmetry 

sin 2 9 

A tt = \S Tp \ |S TS | 1 - cQg2 Q cos(2(/. - <j> Sp - <j> Sp )c4(x p ,Xp) , (15) 

by inserting Eq. (|14fl into Eq. JSJ). In fact, in that case the measurement of A TT would allow the extraction of 
unambigous information on the analytical form of the tranversity h±(x). 

The azimuthal asymmetry defined in Eq. I|5|l. i.e. by flipping the transverse polarization of one of the two beams, 
can be obtained by changing the sign of the cosine function in Eq. (|15fl . While in the laboratory frame the azimuthal 
angles of the beam transverse polarization are fixed, in the Collins-Soper frame they are variable, since the h axis 
is directed along g T /|q T |. Hence, for each randomly distributed </> s _ two sets of events are accumulated for each bin 
in Xp, corresponding to 4> s _ = <f> Sp (parallel transverse polarizations, positive cosine function indicated by U) and to 
(j> s — 4> Sp +tt (antiparallel transverse polarizations, negative cosine function indicated by D). Then, the asymmetry 
is constructed as (U — D)/(U + D) and binned in x p after integrating upon q T , Xp, and the zenithal angle 9. As in 
Ref. the 9 angular distribution is restricted to the range 60-120 deg, because the sin 2 9 dependence would dilute 
the spin asymmetry if events at small 9 were included. This cutoff produces a reduction of events by a factor w 2/5: 
larger statistics are obtained at the price of smaller absolute sizes of the resulting asymmetry. For small modifications 
of the above range, the relative size of statistical error bars and asymmetries will not change much. 

A further reduction of the event sample is due to the transverse polarization of the (anti)proton beams, which is 
assumed to be 50% on the average, giving an overall dilution factor 0.25. In our simulation, this means that an average 
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75% of events have been sorted assuming no polarization at all, while for 25% of events a full transverse polarization 
is assumed for both beams. We hope that actual polarizations will be larger than 50%, but at the same time we must 
be aware that this fact could be compensated by more realistic parton distributions that are less close to the Soffer 
bound than the test functions discussed in this paper. 

III. RESULTS 

In this section, we present results for the Monte Carlo simulation of double spin asymmetries for the p* p^ — » /i + /i~ X 
process in order to explore under which conditions the transversity could be unambiguously extracted from such data. 
As explained in the previous section, the most convenient kinematical option for the HESR at GSI seems, at present, 
the collision of an antiproton beam with energy Ep — 15 GeV and of a proton beam with E p = 3.3 GeV, such 
that the available cm square energy is s w 200 GeV 2 . But we have considered also the option of an antiproton 
beam with the same energy hitting a fixed proton target such that s « 30 GeV 2 . The overall dilution due to beam 
transverse polarization is assumed 0.25. Events are sorted according to the cross section of Eq. JjJJ and supplemented 
by Eqs. |TT |) -l(H |) with \S Tp \ = \S Tp \ = 1 for the 25% of events, and for the 75% of them. 

In the following, we will study the double spin asymmetry of Eq. i|15[l generated by the cos(2(/>— <j) s _ — <fi s ) dependence 
of Eq. (|12|l in the Collins-Soper frame: positive values of the cosine function (U) correspond to parallel transverse 
polarizations of the two beams, negative values (D) to antiparallel polarizations. The asymmetry (U — D)/(U + D) 
is then constructed for each bin x p by integrating upon the other variables q T , x p , 8. In some cases we also show the 
corresponding bidimensional (x p ,x p ) distributions. 

Each displayed histogram contains 17000 events. According to the reduction factors due to the cuts in M, q T , and 
9 (as discussed in the previous section), we need to consider an initial sample of 80000 events, which become 40000 
by applying the cutoff in q T , and are further reduced to the final 17000 after the cutoff in 9. The discussed slightly 
smoother cut q T > 0.5 GeV/c implies a starting sample of 68000 events to arrive at the same final 17000 events of the 
histogram. The q T distribution of Eq. I|1U|I is phenomenological, therefore these reduction factors should be realistic 
(at least in the mass range 4 < M < 9 GeV). In our Monte Carlo simulation it is not possible to predict how many 
muon pairs would be produced with no mass cuts at all, since we take into account physical devices that dominate in 
certain mass ranges only. 

The asymmetry (U — D)/(U + D) will be calculated only for those x p bins that are statistically significant, namely 
that contain a minimum number of events in order to avoid large fluctuations of purely statistical origin. For the 
bidimensional (x p ,x p ) distribution, the cutoff is of at least 10 events per bin. For the Xp-integrated distribution, 
the cutoff is of at least 100 events per bin. Anyway, these statistical cuts do not affect much the overall number of 
surviving events. This latter sample will be referred to as the sample of "good" events, in the sense that it contains 
all the events surviving all the described cuts. As already anticipated, in the histograms of the following figures the 
"good" events amount to 17000. In Ref. |35j, we already discussed the relation between the number of "good" events 
and the running time for an experiment at a given machine luminosity. The same conclusions apply here, so that 
a sample of 10000-30000 events seems a reasonable estimate for a sensible measurement, as it will be clear in the 
following. 

As in the previous paper |35j | , statistical errors are obtained by making 20 independent repetitions of the simulation 
for each considered case, and then calculating for each x p bin the average value of the double spin asymmetry and its 
variance. Again, we checked that 20 repetitions are a reasonable threshold to have stable numbers, since the results 
do not change significantly when increasing the repetitions from 6 to 20. 

In Fig. the sample of 17000 "good" events for the p^ p^ — > /i + /i~ X process is displayed in the above kinematic 
conditions for the collider mode. The invariant mass of the lepton pair is constrained in the range 4 < M < 9 
GeV. The four panels correspond to the choices: a) (hi(x p )} / (fi(x p )) — ^fx^; b) (hi(x p )) / (fi{x p )) = ^/l — x p ; c) 
(hi(x p )} I {fi(x p )) = 1; d) (hi(x p )} / (fi(x p )) = 0. The first two cases have been selected in order to have two opposite 
behaviours, namely ascending and descending, and to verify if they can be identified also in the corresponding 
asymmetries. The last two ones are displayed to set a reference scale for the absolute size of the asymmetry and to 
crosscheck the statistics when A TT cx (hi(x p )) / (fi(x p )} = 0. In any case, all four choices respect the Soffer bound [50| . 
For each bin, two groups of events are stored corresponding to positive values of cos(20 — <fr s — </> Sp ) in Eq. l(T2*|l . 
represented by the darker histograms (events U), and to negative values, corresponding to the superimposed lighter 
ones (events D). The bins at the boundaries, corresponding to x p < 0.1 and x p > 0.7, contain a number of events 
below the discussed statistical cutoffs and they will be discarded in the corresponding asymmetry. 

In Fig. the corresponding double spin asymmetry (U — D) /(U + D) is displayed again in x p bins for all the 
four choices: full squares for (hi(x p )} / (fi{x p )) — y/x^, upward triangles for (h\{x p )) / (fi{x p )) = y'l — x p , downward 
triangles for (hi(x p )} / (fi{x p )) = 1, and open squares for (hi(x p )) / (fi{x p )} = 0. The error bars represent statistical 
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errors only. 
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FIG. 3: The sample of 17000 events for the p T p T — > /i + fj,~ X process where a transversely polarized antiproton beam with 
energy Ep — 15 GeV collides on a transversely polarized proton beam with E p = 3.3 GeV producing muon pairs of invariant 
mass 4 < M < 9 GeV (for further details on the cutoffs, see text), a) (hi(x p )) / (fi(x p )) = ^/x~p~ (brackets mean that each flavor 
contribution in the numerator is replaced by a common average term, similarly in the denominator; for further details, see 
text), b) (hi(xp)) / (fi(xp)} = y/1 — x v . c) (hi{x p )) / (fi(x p )) = 1. d) (h\{x p )) / (fi(x p )) = 0. For each bin, the darker histogram 
corresponds to positive values of cos(2</> — (f> s _ — (f> s ) in Eq. (1121 . the superimposed lighter one to negative values. 
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FIG. 4: Asymmetry (U — D)/(U + D) between cross sections in the previous figure corresponding to darker histograms (U) and 
superimposed lighter histograms (D), as bins in x p . Full squares for the case when {hi(x p )) / (fi(x p )) = ^/x^, upward triangles 
when it equals ^/l — x p , downward triangles when it equals 1 and open squares when it equals 0. Continuous lines are drawn 
to guide the eye. Error bars due to statistical errors only, obtained by 20 independent repetitions of the simulation (see text 
for further details). 

From the above results, we first deduce that the considered sample of events in the specified kinematics is sufficient 
to produce an average significant asymmetry, at most about 15% for (h\{x p )) / (fi(x p )) = 1. When the ratio equals 0, 
the corresponding asymmetry consistently oscillates around in all considered bins, also for low x p where the smaller 
error bars are due to the more dense population of events. In fact, the cross section is known to rapidly increase for 
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decreasing r and data accumulate in the part of the phase space corresponding to the lowest possible r which, for the 
considered kinematics, is r > 0.08. At the same time, in the range 0.1 < x p < 0.4 it seems also that the asymmetries 
corresponding to the ascending ^Jx^ and descending ^/l — x p functions keep this different behaviour. This means 
that, although in this limited (but significant) range, it should be possible to extract also information on the analytical 
dependence of the transversity upon the parton fractional momentum. For higher x v (and r < 0.4), the phase space 
is less populated and the error bars are so large that each one of the four choices for the ratio (hi(x p )) / (fi(x p )) can 
be confused with the other ones. 




FIG. 5: Unintegrated asymmetry (U — D)/(U + D) for the case (hi(x)} / (fi(x)} = *Jx in the same conditions as the previous 
figure, but plotted in bins of Xp and x p . Left panel: distribution of average values. Right panel: distribution of the variances, 
i.e. of half the statistical "error bars", obtained by 20 independent repetitions of the simulation (see text for further details). 

This trend is confirmed and better clarified by looking at the unintegrated asymmetry, displayed in Fig. [S] in bins 
of Xp and x p for the case (hi(xp)) (hi(x p )} / (fi{xp)) (fi(x p )} = ^/x^ y/x~p~. Namely, it corresponds to the full squares 
of Fig. for the ^-integrated case. The left panel represents the tridimensional plot of the average values of the 
unintegrated double spin asymmetry, while the right panel gives the distribution of the variance, i.e. of half the 
"error bar" for each (x p ,x p ) bin. On the boundaries of the plot, where close to 1, the small number of 

collected events produces large statistical errors and also large asymmetries due to fluctuations. Anyway, each bin 
containing less than 10 events is considered not statistically relevant and the corresponding asymmetry and error have 
been artificially put to zero. Viceversa, for small x p < 0.3 the variance is small through all the x p range so to give a 
distinguishable integrated double spin asymmetry. Unfortunately, because of the ascending trend of the function a/x, 
the absolute values of the asymmetry are small in the x p range of interest, even if statistically different from zero. 




FIG. 6: Unintegrated asymmetry (U — D)/(U + D) for the case (h\(x)) / (fi{x)) = VI — x in the same conditions as the previous 
figure, plotted in bins of Xp and x p . Left panel: distribution of average values. Right panel: distribution of the variances, i.e. 
of half the statistical "error bars", obtained by 20 independent repetitions of the simulation (see text for further details). 

In Fig. 0the unintegrated double spin asymmetry is displayed for the case (hi(xp)) (hi(x p )) / (fi(xp)) (fi(x p )) = 
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yj\ — Xp — x p in the same conditions as in the previous figure. Therefore, it corresponds to the upward triangles 
of Fig. 01 for the Xp-integrated case. Similar arguments apply to the statistical selection of the results. The only 
difference is that for the relevant x p < 0.3 range the absolute values of the asymmetry are more significant because of 
the descending trend of the function yl — x. By comparing the right panel of this figure with the corresponding one 
in Fig. [SJ we notice that errors actually do not depend on the selected test function. They increase with increasing r, 
with the exception of some bins that are crossed by the hyperbole M > 4 GeV/c. 



events 




FIG. 7: The sample of 17000 events for the p T p T — > /i + fj,~ X process for the lepton invariant mass 1.5 < M < 2.5 GeV, while 
the other kinematic conditions, conventions for each panel, and color codes of histograms are as in Fig|3] (for further details, 
see text). 

In Fig. Q, the sample of 17000 "good" events for the — > /i + fj," X process is displayed in the same conditions 
as in Fig. |3 but with the invariant mass of the lepton pair in the range 1.5 < M < 2.5 GeV. Again, the conventions for 
the four panels and the color codes of the histograms are the same as in Fig. 13 Since the range of explored r = M 2 /s 
is now 0.01-0.03, the bins for x p — > are more populated than before, while for x p > 0.7 yet they are statistically not 
significant and again they will be discarded in the asymmetry plot. 

In Fig. |SJ the corresponding double spin asymmetry (U — D)/(U + D) is displayed with the same conventions 
as in Fig. 21 The larger population for small x p bins reflects in smaller error bars for all four choices of the ratio 
(h\{xp)) / (fi(x p )) . Consequently, in the range 0.1 < x p < 0.3 the functions ^~x~p~ (full squares) and yjl — x p (upward 
triangles) are even more distinguishable than in the previous case. 

As already recalled in the previous section, we have also explored the typical kinematics suggested by the PANDA 
collaboration in its proposal at HESR at GSI [4(|. Namely, the operational mode where the antiproton beam with 
energy E p = 15 GeV hits a proton target and produces lepton pairs with low invariant masses; we have considered 
the range 1.5 < M < 2.5 GeV for the same reasons as above. We have still assumed that the transverse polarization 
of both beam and target is 50% on the average, such that the corresponding overall dilution factor is 0.25. This 
corresponds to take |S T | = |S r J = 1 in Eqs. 1)12(1 and 1)15(1 for the 25% of the events, while taking them for the 
remaining 75%. The resulting cm square energy is s = 30 GeV 2 and the range of explored r is approximately 0.07-0.2. 
Consequently, the integrated event distribution of Figgis more populated at higher x p bins for all the four 
ratios explored (notations, conventions and histogram color codes are as in Fig. 01 . 

In Fig. ^3 the double spin asymmetry (U — D)/(U + D), corresponding to cross sections of the previous figure, 
is displayed with the same conventions as in Fig. 21 At variance with the result of Fig. |H| the explored portion of 
phase space for larger r (hence, for larger x p bins) is not favoured by the given 1/r qualitative dependence of the 
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FIG. 8: Double spin asymmetry (U — D)/(U + D) between cross sections in the previous figure corresponding to darker 
histograms (U) and superimposed lighter histograms (D), as bins in x v . Notations as in Fig.|3] Continuous lines are drawn to 
guide the eye. Error bars due to statistical errors only, obtained by 20 independent repetitions of the simulation (see text for 
further details). 
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FIG. 9: The sample of 17000 events for the p T — + fi + /i~ X process where a transversely polarized antiproton beam with 
energy Ep — 15 GeV hits a transversely polarized proton target producing muon pairs with invariant mass 1.5 < M < 2.5 GeV 
and s = 30 GeV 2 . Conventions for each panel and color codes of histograms are as in Fig|^l(for further details, see text). 

cross section: a lower density of events reflects in larger statistical error bars which prevent from clearly distinguishing 
each one of the four considered forms for the (/ii)/(/i) ratio. This means that for the considered sample of 17000 
"good" events, it is not yet clear how to extract information on the analytical structure of hi(x). However, it should 
be possible at least to observe a nonvanishing asymmetry and give an estimate of its magnitude and sign. 

The relevant message of previous figures is that it is crucial to consider integrated distributions in one parton 
fractional momentum only, in order to reasonably populate bins and to reach a deconvolution of transversity from the 
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FIG. 10: Double spin asymmetry (U — D)/(U + D) between cross sections in the previous figure corresponding to darker 
histograms (U) and superimposed lighter histograms (D), as bins in x p . Notations as in Fig. El Continuous lines are drawn to 
guide the eye. Error bars due to statistical errors only, obtained by 20 independent repetitions of the simulation (see text for 
further details). 

product h(xp)h(x p ). The price to pay is that the useful phase space in the other fractional momentum is reduced to 
few bins. It is easy to estimate where the maximum of the distribution is located. Assuming that the bidimcnsional 
(x^, x p ) distribution is dominated by the 1/r factor associated with the elementary qq fusion into a virtual photon, 
the distribution is l/(xpx p ) for r > M% lin /s, and otherwise. Consequently, the ^-integrated distribution has the 
form \og(x p s/M? nin )/xp and reaches its peak value for \og{x p s / M^ lin ) = 1, i.e. for x p = eM^ in /s « 3M? nin /s. For 
M = 4 GeV and s — 200 GeV 2 , x p « 0.22 in agreement with Fig. It is evident that it is possible to span the 
domain of valence contributions for s in the range 100-300 GeV 2 and M > 4 GeV. A similar conclusion holds by 
directly considering the unintegrated distribution, where the relevant contribution of annihilating valence (anti)partons 
corresponds to Xp — x p w 0.3, such that t — x p xp ss 0.1 can be reached again for 100 < s < 300 GeV 2 and M > 4 
GeV. 

Moreover, exploring different masses at the same s (as we did in Figs. 0] and |8J) is useful to estimate the role of 
higher twist effects, since the latter can be classified according to powers of M/M p , where M p is the proton mass. 
However, our results do not include these corrections since precision calculations are beyond the scope of the present 
work. 

In conclusion, it seems that in the collider mode for the HESR at GSI, where the cm square energy is s = 200 
GeV 2 , a sample of 17000 Drell-Yan events from collisions of transversely polarized antiproton and proton beams with 
proper lepton invariant masses is sufficient to generate sizeable double spin asymmetries from which information can 
be extracted about the functional dependence of the transversity, but limited to the range 0.1 < x < 0.3. For the 
lower case s = 30 GeV 2 in the fixed-target mode, the double spin asymmetry is still sizeable, but the larger statistical 
error bars do not allow for such a clean extraction. 



IV. CONCLUSIONS 

In a previous paper [35| , we produced a Monte Carlo simulation to study the physics case of a Drell-Yan process with 
unpolarizcd antiproton beams and transversely polarized proton targets. Here, we have considered the fully polarized 
pT p\ _^ ^+ jj- x process in order to explore the leading transverse spin structure of the nucleon. Both works are 
finalized at the kinematics of the High Energy Storage Ring (HESR), a source for (polarized) antiprotons under 
development at GSI. In fact, using antiproton beams offers the advantage of involving unsuppressed distributions of 
valence partons: the transversely polarized quark in a transversely polarized proton, and the transversely polarized 
antiquark in a transversely polarized antiproton. 

Different kinematical options have been considered. Since the cross section fastly decreases for increasing r = 
M 2 /s = XpXp (where M is the lepton pair invariant mass and s is the center-of-mass square energy), events tend to 
accumulate in the phase space part corresponding to the smallest r allowed by the mass cutoff, which is dominated 
by the valence contribution to parton distributions. It seems convenient to reach such low values of r by adequately 
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increasing s. The Xp-integrated distributions present a Poisson-likc shape in the parton fractional momentum x p with 
a peak at x p w 3M^ in /s and width 0.1-0.2. Outside this range, the event distributions arc of difficult analysis. For 
typical M min values, the most interesting s range turns out to be 100-300 GeV 2 , where the possibility of repeating 
the experiment at different s values allows to analyze different x p ranges. 

We have simulated the fully polarized Drell-Yan process using antiproton beams with energy E p = 15 GeV and 
(transverse) polarization 50%, and proton beams with the same polarization and energy E p = 3.3 GeV, such that 
s = 200 GeV 2 . To avoid the cc threshold and other resonances in the mass spectrum of the lepton pair, we selected 
the invariant mass M in the ranges 4-9 GeV and 1.5-2.5 GeV. The corresponding explored r ranges are 0.08-0.4 
and 0.01-0.03, respectively. We have also explored the case where the transversely polarized antiproton beam of 15 
GeV hits a transversely polarized fixed proton target producing lepton pairs with invariant mass in the same range 
1.5-2.5 GeV. In this case, s = 30 GeV 2 and the r range is approximately 0.07-0.2. The mass range 4-9 GeV has not 
been considered here, because it implies average x p values where the parton model underlying the simulation is not 
appropriate. 

In all cases, the transverse momentum of the dimuon couple has been limited to q T > 1 GeV/c and its zenithal- angle 
distribution is restricted to the range 60-120 deg. The former cut is induced by the need to avoid complicated soft 
mechanisms and to match the experimental requirements of the collider setup. The latter cut prevents the angular 
dependence of the cross section from diluting the related asymmetry. All these cuts produce a remarkable reduction 
of the Drell-Yan events: the considered initial sample of 80000 events (only affected by the mass cut) is reduced to 
17000. 

At leading twist in the cross section, the contribution of interest has a characteristic azimuthal dependence of the 
kind cos(2</> — <j) s _ — <fi Sp ), where <ft is the azimuthal angle of the final lepton pair and <p s is the azimuthal position 

of the (anti)proton transverse spin in the Collins-Soper frame. Hence, for each randomly distributed </> s _ events have 
been accumulated for parallel {4> s _ = (f> Sp ) and antiparallel ((f> s _ — <p Sp + w) beam polarizations. The corresponding 
azimuthal asymmetry has been built either integrating upon all variables but the quark fractional momentum x p , or 
keeping also the corresponding antiquark x p dependence and building bidimcnsional plots. 

For each kinematical case, we have considered four functional dependences for the transversity, each one with a 
very peculiar trend, namely constant, ascending, descending, and zero, but all satisfying the Soffer bound. The goal 
is to recover the same different trend also in the corresponding double spin asymmetry, which means that feasibility 
conditions for an unambiguous extraction of the transversity could be established. In the collider mode with s = 200 
GeV 2 , the double spin asymmetry has very small statistical error bars for 0.1 < x p < 0.3. With this limitation, it 
seems that a sample of 17000 events satisfying all the described cutoffs, should be sufficient to grant the identification 
of the analytical behaviour of the transversity. This statement is still valid at the lower M range considered by 
keeping the same s, because the resulting r range is shifted to lower values and the statistics is higher. Viceversa, for 
the option where the antiproton beam hits a fixed proton target, a lower s results in larger error bars: measuring a 
nonvanishing double spin asymmetry seems possible, but the extraction of the transversity looks more problematic. 

In conclusion, with the present simulation we have explored the feasibility conditions for an unambiguous extraction 
of the transversity from Drell-Yan data with polarized antiproton beams; we hope to have contributed to the studies 
of the physics case for hadronic collisions with antiproton beams at the HESR at GSI. 
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